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— ‘  For  an  arbitrary  conducting  surface  Z(jc,  y)  we  obtained  an  analytical  expression  for  the  local  refractive  index >  as  a 
function  of  Z,  dZ/dx,  dZ/Qy,  and  the  Drude  conductivity  *  by  using  the  complex  ray-tracing  method.  The  Fresnel 
coefficients  of  reflectance  and  transmittance  are  then  employed,  and  the  value  of  y  is  obtained  to  determine  the 
scattered  and  refracted  fields.  The  proposed  method  has  advantages  over  the  methods  of  solution  by  the  Debye 
potential,  the  Laplace  transform,  and  the  vector-wave  equation  in  the  computation  of  the  scattering  and  absorption 
parameters  of  a  wide  range  of  complex  surface  and  wave-front  geometries.  The  surface  integrals  obtained  in  the 
present  study  include  the  surface  function  in  an  explicit  and  concise  form.  t  ■  , 


1.  INTRODUCTION 

Since  Mie  (1908)  introduced  his  rigorous  solution  of  Max¬ 
well’s  equations  for  the  diffraction  of  a  plane  monochromat¬ 
ic  wave  by  a  conducting  sphere  by  using  Debye’s  potentials,  a 
great  number  of  applications  have  been  made  by  using  his 
approach  in  solving  diffraction  problems  involving  particles 
of  various  shapes.1  The  main  drawback  of  using  the  Debye 
potential  method  is  the  difficulty  of  finding  suitable  solu¬ 
tions  for  the  scalar-wave  equation  to  represent  either  an 
incident  wave  front  other  than  a  plane  wave  or  the  scattered 
and  refracted  waves  in  complex  geometries.  Another  ap¬ 
proach  to  the  solution  of  the  diffraction  problems  is  to  solve 
the  vector-wave  equation  in  the  space  external  to  the  bound¬ 
ary  surface.  The  solution  obtained  in  this  way  is  a  function 
of  the  integral  of  the  spherical  vector  harmonics  at  the  sur¬ 
face.2  A  third  approach,  in  the  long-wavelength  range,  is  the 
solution  of  the  Stratton-Chu-Silver  integral  (SCSI)  form 
with  the  perfectly  conducting  surface  and  tangent  plane 
approximations.3-4 

In  the  visible  range  of  wavelength  the  conductivity  is  fi¬ 
nite.  Leader*  worked  out  the  scattering  by  a  rough  surface 
and  obtained  a  solution  of  the  SCSI  form  with  the  deriva¬ 
tions  of  the  Fresnel  coefficients  not  expressed  explicitly  in 
terms  of  the  surface  function  Z(x,  y).  Also,  Bahar  and 
Rajan9  analyzed  the  same  problem  by  using  a  complicated 
mathematical  formulation. 

In  the  present  paper  we  present  a  solution  palatable  to 
workers  in  the  field  of  optics.  We  have  utilized  the  general¬ 
ized  ray-tracing  procedure  developed  by  Kneisly7  and  oth¬ 
ers8'9  to  analyze  the  diffraction  of  the  electromagnetic  waves 
at  an  arbitrary  surface  of  a  finite  conductivity.  The  defini¬ 
tion  of  the  refractive  index  for  a  conducting  medium  used  by 
Stratton10  for  a  plane  surface  is  generalized  for  an  arbitrary 
surface.  The  Fresnel  coefficients  of  reflectance  and  trans¬ 
mittance  are  then  employed  with  the  new  geometry-  and 
medium-dependent  refractive  index  in  the  place  of  the  tra¬ 
ditional  only-medium-dependent  refractive  index. 

We  will  ignore  internal  reflection  in  the  conducting  medi¬ 
um,  since  scattering  is  greatly  overwhelmed  by  absorption. 

_ This  simplifies  the  boundary  conditions  to  a  great  extent. 
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Our  analysis  is  also  restricted  to  geometries  in  which  no 
multiple  reflection  occurs  between  opposing  surfaces.  As 
we  will  see  in  the  application  of  the  present  method,  the  Mie 
problem  is  a  good  candidate  for  our  method  of  solution 
because,  for  a  metallic  sphere,  it  is  realistic  to  neglect  br  *h 
internal  and  external  multiple  scattering.  Treatment  of 
penetrable  and  inhomogeneous  bodies  was  carried  out  by 
others  for  simpler  surface  geometries  by  using  different 
methods  of  numerical  computation.11-13  A  study  was  done 
for  refracted  fields  from  an  arbitrary  refractive  surface,14  but 
reflection  was  ignored  in  that  study. 

In  the  applications  in  which  the  value  of  the  field  at  the 
surface  is  the  main  concern,  such  as  in  laser-solid  interac¬ 
tion,  our  solution  is  expressed  in  terms  of  modified  Fresnel 
coefficients.  In  the  applications  in  which  the  scattered  field 
is  desired,  such  as  in  radar  measurements  and  remote  imag¬ 
ing  applications,  the  solution  is  expressed  in  a  modified 
SCSI  form. 

The  rest  of  the  paper  falls  into  six  sections.  In  Section  2, 
we  obtain  an  analytical  expression  of  the  refractive  index  v  in 
terms  of  the  surface  Z(x,  y),  Zx,  Zy,  and  the  physical  proper¬ 
ties  of  the  media.  In  Section  3  we  determine  the  different 
electromagnetic  field  vectors  of  the  locally  polarized  waves 
on  the  surface  Z(x,  y).  In  Section  4  the  absorption  of  the 
incident  energy  is  analyzed.  In  Section  5  the  scattering 
fields  are  described.  In  Section  6  an  application  to  a  spheri¬ 
cal  wave  front  converging  on  a  finitely  conducting  surface  is 
discussed.  Our  conclusion  is  stated  in  Section  7. 

2.  REFRACTIVE  INDEX  OF  AN  ARBITRARY 
CONDUCTING  SURFACE 

The  Stratton  definition  of  the  refractive  index  of  a  conduct¬ 
ing  plane  surface  can  be  generalized  to  a  surface  of  an  arbi¬ 
trary  coordinate  function  as  follows.  At  a  surface  of  discon¬ 
tinuity  Z(x,  y)  separating  two  media,  the  local  normal  to  the 
surface,  fl,  makes  a  polar  angle  6  with  the  z  axis  and  an 
azimuth  with  the  x  axis  in  the  xy  plane  of  the  reference 
frame.  These  angles  are  given  by  differential  geometry 
principles  as  (see  also  Ref.  3  and  4) 

©  1988  Optical  Society  of  America 

88  3  14  1  53 


M.  F.  El-Hewie  and  R.  J.  Cook 


Vol.  5,  No.  2/F ebruary  1988/J.  Opt.  Soc.  Am.  A  201 


tan  0  =  ~(ZX 2  +  Zv2)1/2,  (la) 

tan  y?  =  Zv/Zx,  (lb) 

where 

Zx  =  dZ/dx ,  Zv  =  0Z/dy.  (lc) 

The  subscripts  x  and  y  denote  differentiation  only  when 
they  appear  on  Z;  they  denote  Cartesian  components  else¬ 
where.  The  unit  normal  vector  h, 

h  =  i-Zj  -  Zj  +  k)/C,  (2a) 

is  then  determined,  where 

C=  (1 +ZI2  +  Z>2)1/2.  (2b) 

The  propagation  unit  vector  normal  to  the  incoming  wave 
front  A,  is  given  by 

h,  =  sin  0t  cos  vy  +  sin  6,  sin  s pj  +  cos  dfi,  (2c) 


where  8,  and  are  given  independent  of  the  position  vector 
of  a  point  at  the  surface,  r,  for  a  plane-wave  illumination. 
At  a  point  P(x,  y,  Z )  on  the  surface  Z,  the  unit  vectors  along 
the  refracted  and  reflected  rays  are  of  well-known  form  [see, 
for  example,  Eqs.  (22)  and  (29)  of  Ref.  7]  and  are,  respective¬ 
ly,  given  by 

A,  =  (kjkr)ht  +  [cos(0,„  +  7r)  -  (kjkt) cos  0in]A  (2d) 
and 

A,  =  A,  -  2  cos  0inA.  '  (2e) 

ir  is  added  to  drn  because  the  refracted  rays  are  propagating 
away  from  the  positive  z  direction.  The  subscripts  i,  s.  and  r 
refer  to  quantities  related  to  the  incident,  reflected,  and 
refracted  waves,  respectively.  Double  subscripts  on  angles 
indicate  the  angle  between  the  propagating  vectors  .of  the 
designated  waves  and  the  surface  normal;  e.g.,  6rn  is  the  angle 
between  Ar  and  A.  k,  and  kr  are,  respectively,  the  wave 
propagation  numbers  of  the  incident  and  the  refracted 
waves,  k ,  is  real,  whereas  kr  is  complex.  The  local  angles  of 
incidence,  6 in(x,  y,  Z)  and  <£;n(x,  y,  Z),  are  given  by 

y,  Z)  =  w\„  =  =  v?(x,  y,  Z)  -  <e,(x,  y,  Z),  (3a) 

cos  0in  =  A  •  A , 

=  (-Zj  sin  0,  cos  <e,  -  Zy  sin  0,  sin  +  cos  0;)/C 

(3b) 

and 

sin  fljn  =  (I  -  (A  •  n,)2]1/2 

=  |I  -  ((— Zr  sin  0,  cos  v,  -  Z,  sin  0,  sin 
+  cos0,)/C]2|l/2.  (3c) 

For  the  refracted  waves,  Snell’s  law  in  its  scalar  form  gives 
the  local  angle  of  refraction,  0r„(x,  y,  Z),  as 

kr  sin  0rn  =  A,  sin  0in.  (4) 

The  arguments  x,  y,  and  Z  are  sometimes  suppressed  to 
simplify  notation  whenever  there  is  no  chance  for  confusion. 
An  aplanar  wave  front  is  described  by  a  ray-aberration 


vector  u(ux,  uy,  u,).  The  components  of  u  are  the  projec¬ 
tions  of  the  phase  shift  along  an  incident  ray  on  the  different 
Cartesian  directions.  Now,  consider  the  following  general 
forms  for  the  incident,  refracted,  and  scattered  waves1; 

e,(r,  t)  =  e;(r)  exp[j(fe(A,  -  r  +  A,u  •  r  -  u>£)J,  (5a) 

h,(r,  t)  =  h,(r)  expfRA.A,  •  r  +  k,u  •  r  -  w£)],  (5b) 

er(r,  t)  =  e,(r)  exp[i(A,.A,.  •  r  +  kp  •  r  -  «£)],  (5c) 

hr(r,  £)  =  hr(r)  exp[i (k/ir  •  r  +  A,u  •  r  -  «t)],  (5d) 

e,(r,  £)  =  e,(r)  exp[t(A,A,  •  r  +  A,u  •  r  -  &>£)],  (5e) 

h,(r,  t)  =  h,(r)  exp(((A,A,  •  r  +  A,u  •  r  -  <at)],  (5f) 

with  e(r)  and  h(r)  terms  real.  er(r,  £),  hr(r,  t),  e,(r,  t),  and 
h,(r,  t)  are  the  refracted  and  reflected  fields  and  are  related 
by  the  Maxwell  curl  equations, 

Mr*1  r(r>  t)  =*  -V  X  er(r,  t )  (6a) 

and 

M,[J»,(r,  t)  +  h,(r,  t))  =  -V  X  [e,(r,  t)  +  e,(r,  t)).  (6b) 

On  using  Eqs.  (5)  in  Eqs.  (6),  with  some  manipulation,  we  get 
the  following  two  expressions  for  the  refracted  and  reflected 
magnetic  fields: 

hr  *  V(k,hr  •  r/ki  +  u  •  r)  X  ejv,  (7a) 

h,  *  h,  +  [V(A,  •  r  +  u  •  r)  X  e,  +  V(A,  •  r  +  u  •  r)  X  e,]/q 

(7b) 

since  the  curl  of  the  amplitudes  of  the  incident,  reflected, 
and  scattered  fields  vanish  at  the  boundary  surface.  The 
angular  frequency  of  light  u  *  and  the  intrinsic 

impedance  of  free  space  ij  =  Oq/q)*''2  are  used.  m.  and  t,  are, 
respectively,  the  magnetic  permeability  and  the  electric  per¬ 
mittivity  of  the  free  space. 

The  real  and  imaginary  parts  of  the  refracted  wave  propa¬ 
gation  vector  krhr  are  obtained  as  follows.10  kr  and  cos  0rn 
are  written  explicitly  ir  •  heir  complex  forms  as 

kr  *  £  +  i{  (8a) 


cos  0rn  *  Sk  +  iA,  (8b) 

where  £  and  f  are  known  from  the  elementary  electron  the¬ 
ory  of  the  optics  of  metals  and  are  conveniently  given  by1 

£  *=  <40.5„r<r((l  +  (jr2/(r2u)2)1/2  +  1]|1/2  (8c) 


f  =  w!0.5^r<r[(l  +  -  1]|1/2.  (8d) 

or  is  the  Drude  electrical  conductivity  of  the  metal.  ♦  and  A 
are  functions  of  Z,  Z„  Zy,  0„  and  the  physical  properties  of 
the  media  and  are  obtained  from  Eq.  (4).  After  squaring  Eq. 
(4),  using  Eq.  (8a),  and  performing  some  manipulation,  we 
get 


*  *  |0.5[(A2  +  B2)1'2  +  A]|1/2 


*vtui  axi4/«r 
|  Spaolal 
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A  =  |0.5[(A2  +  B2)'n  ~  A]]m,  (80 

where 

A  =  1  -  k2  sin2  0,n({2  -  ^)/(£2  +  f2)2  (8g) 

and 

B  =  2*l2(ff)sin2  6j(k2  +  ?)2.  (8h) 

After  substituting  sin  0-,n  from  Eq.  (3c)  into  Eqs.  (8g)  and 
(8h),  we  get 

A  =  1  -  A,2[($2  -  ^)/((2  +  ^)2]|1  -  [(-2,  sin  6 ,  cos 

-  Zy  sin  0,  sin  +  cos  0,)/C]2|  (8i) 

and 

B  =  2A;2[£f/(£2  +  f2)2]!!  -  l(-Zx  sin  0f  cos  v?4 

—  Zy  sin  0,  sin  w,-  +  cos  0;)/C2|.  (8j) 

After  substituting  cos  0rn  from  Eq.  (8b)  into  Eq.  (2c),  we  get 

hr  =  (kjkr)hl  -  [*  +  i A  +  (kjkr)  cos  0in]A.  (9a) 

By  multiplying  both  sides  of  Eq.  (9a)  by  kr  from  Eq.  (8a)  and 
then  arranging,  we  get 

Mr  =  *A  “  (W  “  fA  +  *,•  cos  0in)A  -  i(f*  +  ?A)A 

=  ^,[n,  +  (q  +  ip)A].  (9b) 

After  substituting  cos  dm  from  Eq.  (3b)  into  Eq.  (9b),  we  get 
q  =  -({'I'  -  fA)/fe,  +  (Zx  sin  0;  cos  ^  +  Zy  sin  6t  sin 

—  cos  0j)/C  (10a) 

and 

p  «  -({*  +  {A)/fcj.  (10b) 

The  q  and  p  defined  by  Eqs.  (10)  for  an  arbitrary  surface  are 
similar  to  those  obtained  by  Stratton  for  a  plane  surface,  but 
they  are  now  functions  of  2  and  its  derivatives.  After  sub¬ 
stituting  Mr  from  Eq.  (9b)  into  Eq.  (7a),  we  get 

h,  -  V[(A,  +  (q  +  ip)h  +  u)  •  r]  X  e^ij.  (11) 

Equation  (11)  relates  the  refracted  magnetic  and  electric 
fields  algebraically,  since  the  gradient  term  is  only  a  function 
of  the  surface  geometry  and  parameters.  For  an  incident 
plane  wave,  A,  depends  on  r  through  h  [Eq.  (2d)].  For  an 
arbitrary  wave  front,  Ar  depends  on  r  through  both  A,  and  A. 
Equation  (11)  is  a  direct  result  of  our  assumption  that  only  a 
single  refracted  wave  exists  inside  the  conducting  medium, 
which  simplifies  the  boundary  conditions  and  makes  it  pos¬ 
sible  to  determine  the  surface  values  of  the  field  algebraical¬ 
ly.  Our  assumption  is  justified  on  the  basis  that  the  skin 
depth  in  metals  is  ultimately  small  compared  with  the  wave¬ 
length  of  the  field.  This  is  not  the  case  with  dielectrics, 
where  e,  is  a  sum  of  an  ingoing  wave  and  an  internally 
reflected  wave. 

Equation  (11)  represents  a  system  of  inhomogeneous 
waves.10  The  imaginary  part  of  the  gradient  term  in  Eq. 
(11)  represents  the  attenuation  of  the  conducting  medium. 
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The  real  part  of  the  gradient  term  is  the  local  wave  propaga¬ 
tion  vector  along  the  refracted  ray  and  can  be  written  as 

K  =  A,  +  qh  +  u.  (12) 

The  real  local  angle  of  refraction  0,/  is  then  defined  as 

cos  0r(  =  K  •  A/lKl  (13a) 

and 

sin  0,(  =  [1  —  (K  •  A/|K|)2]I/2.  (13b) 

By  substituting  A,  and  A  from  Eqs.  (2)  and  K  from  Eq.  (12) 
into  Eqs.  (13),  we  get 

cos  erl  =  {-ZXF,  -  ZyF2  +  F3  +  qC)/[(CF ,  -  qZx)2 

+  (CF2-qZy)2  +  (CF3  +  q)2]U2  (13c) 


and 


sin  8ri 

-  [1  -  (-ZXF,  -  ZyF2  +  F3  +  qC)2/[(CFl 

-  qZx )2 

+  (CF2  -  qZy)2  +  (CF3  +  q)2]]1/2, 

(13d 

where 

* 

F,  =  sin  0,  cos  v>,  +  ux, 

*  (14a) 

and  • 

F2  ■  sin  0,  sin  +  uy, 

■  (14b) 

F3  *  cos  0j  +  ux. 

(14c) 

The  u’s  vanish  for  a  plane  wave  and  should  be  given  for  an 
aberrated  wave.  From  Eqs.  (3c),  (13d,)  and  (14),  the  local 
refractive  index  is  given  as 

i/(2,  a,  u)  ■  sin  0in/sin  0r/ 

-[(1  -  |[-ZX(F,  -  ux)  -  Zy(F2  ~  Uy)  +  F3 

-  uJ/C|2)/|l  -  (~ZxFl  -  ZyF2  +  F3 

+  qC)2/[(CF i  -  qZx)2  +  ( CF2  -  qZy)2 

+  (CF3  +  <?)2]|]1/2.  (15) 

Equation  (15)  gives  the  local  refractive  index  in  terms  of  2 
and  its  derivatives,  the  angles  of  incidence  with  respect  to  a 
reference  frame,  the  optical  abjiie ration  u,  and  the  physical 
properties  of  the  media.  Equation  (15)  can  be  reduced  to 
the  case  of  a  plane  wave  irradiating  a  plane-conducting  sur¬ 
face  as  follows.  Let  the  surface  derivatives  Z»  and  Zy  vanish, 
and  let  <fii  m  0  and  u  »  0;  we  get  from  the  above  equations 

F,  ■  sin  0„  F2  “0,  F3  -  cos  0„  C  ■  1, 

q  -  -({¥  -  f A)/*,  -  cos  Bj.  (16a) 

Substituting  the  different  values  from  Eq.  (16a)  into  Eq. 
(15)  gives 

»su.t«o.  “  I*,2  *in2  0,  +  U*  -  fA)2]1'2/*,,  (16b) 

which  was  obtained  by  Stratton10  with  the  term  (f'F  -  fA) 
equivalent  to  his  q. 
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3.  ELECTROMAGNETIC  FIELD  VECTORS  OF 
THE  LOCALLY  POLARIZED  WAVES  ON  THE 
SURFACE  Z[x,  y) 

Lei  a  be  the  polarization  unit  vector  and  let  E,  be  the  inci¬ 
dent  electric  field  strength.  Also  let  a  make  an  angle  a  with 
the  plane  that  contains  both  the  unit  vectors  ft,  (along  the 
incident  ray)  and  £  (along  the  z  axis).  From  an  orthogonal 
system  of  unit  vectors,  a  can  be  written  in  the  following 
Cartesian  form  (see  Fig.  1  and  Refs.  5  and  15) 

a  =  (-cos  a  cos  9,  cos  if,  -  sin  a  sin  if,)i 


+  (-cos  a  cos  9 ,  sin  if,  +  sin  a  cos  if,)j 


+  cos  a  sin  9,ft. 


At  any  local  point  at  the  surface  Z  the  incident  fields  can  be 
decomposed  into  both  a  locally  polarized  s  wave  (TE)  and  a 
p  wave  (TM)  as  follows.  A  local  orthogonal  system  of  unit 
vectors  is  defined  as 


t  =  ( h ,  X  n)/[l  -  (A,  •  h)2 


’  p=  h  x  t, 


t}=  h,  X  t. 


After  substituting  fi,  and  rt  from  Eq.  (2)  into  Eqs.  (17b)— 
(17d)  and  executing  the  vectorial  multiplication,  we  finally 
get 


f  =  ( 1/C  sin  9m)[i(sin  9,  sin  if,  +  Zv  cos  9,) 
+  /(-sin  8 ,  cos  v*,  -  Z,  cos  6,) 

+  £  sin  8,(ZX  sin  if,  -  Zv  cos  v?,)]. 


p  =  ( 1/C2  sin  9,n)|r[sin  9,(Z2  cos  <f,  -  ZxZy  sin  if, 


+  cos  ,z,)  +  Zx  cos  9,j  +  /[sin  9,(Z  2  sin  f, 

-  ZxZy  cos  +  sin  v,)  +  Zy  cos  9,j  +  fc[cos  9,(Z2 


+  Z2)  +  sin  9,(Z,  cos  f,  +  Z ,  sin  <p,)]|, 


a  ■  t  =  (  -1/C  sin  9in)[Z,(sin  a  cos  9,  cos  f, 


-  cos  o  sin  if,)  +  Z,(sin  a  cos  9,  sin  if, 
■f  cos  a  cos  sf,)  +  sin  a  sin  8,], 


a  ■  n  =(l/C)[ZI(cos  a  cos  9,  cos  tf,  +  sin  a  sin  if,) 


+  Z>( cos  a  cos  9,  sin  f,  -  sin  «  cos  y:,) 
+  cos  a  sin  9,], 


d  ■  p  =  ( 1/C2  sin  9jn)|(-cos  a  cos  9,  cos  f, 


-  sin  a  sin  <£,)[sin  9,(ZV2  cos  f,  -  ZxZy  sin  if, 


+  cos  f,)  +  Zx  cos  9,j  +  (-cos  a  cos  9,  sin  if, 

+  sin  a  cos  sc,)  [sin  9,(ZX2  sin  if,  -  ZxZy  cos  if, 

+  sin  f,)  +  Zy  cos  9,)  +  cos  a  sin  9, (cos  9,(Z,2  +  Z2) 
+  sin  9,(Z,  cos  v>,  +  Zy  sin  <c,))|.  (18e) 
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Fig.  1.  Geometry  of  the  irradiated  surface. 


For  an  s-polarized  (TE)  wave,  the  field  vector  components  at 
the  surface  are  then  given  by 


ei  -  (A  •  t)'tE, 


hi  =  ft,  X  ei/rj, 


and  for  a  p-polarized  (TM)  wave  the  field  vector  compo¬ 
nents  at  the  surface  are  given  by 


e‘p  =  [(«  •  h)h  +  W  •  (>)P\E, 


hp  -  h,  X  e'p/rt.  (19d) 

On  e  and  h  the  subscript  s  denotes  scattering,  but  on  e‘  and 
h‘  it  denotes  an  s  wave  (TE).  The  superscripts  i  and  a 
(below)  stand,  respectively,  for  incidence  and  absorption. 
The  amount  of  local  incident  energy  flow  per  unit  area  nor¬ 
mal  to  the  surface  in  each  wave  is  then  defined  as 


</i  =  0.5e‘  •  e*‘  cos  8,Ji\ 


jp  m  0.5ej,  •  e*'  cos  9in/ij. 


By  substituting  cos  9jn  from  Eq.  (3b)  and  the  e'  terms  from 
Eqs.  (19)  into  Eqs.  (20),  we  obtain  expressions  for  the  inci¬ 
dent  power  intensity  in  terms  of  Z  and  its  derivatives  and  the 
angles  of  incidence: 

«/j  «  0.5(E,2  cos  8JVM  ■  t)2 

-  0.5[E,2  cos  9in/»)(C  sin  9in)2] 


X  |Z,(sin  a  cos  8,  cos  if,  -  cos  a  sin  y>,) 

+  ZJsin  a  cos  9,  sin  +  cos  a  cos  if,)  +  sin  a  sin  9,|2, 


^  *■  -Y‘ 


‘ 
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Ji  =  0.5(£,2  cos  0in/ij)[(a  •  A)2  +  (5  •  p)2}. 


(22) 


All  terms  in  Eq.  (22)  were  obtained  above  in  terms  of  the 
angles  of  incidence  and  the  surface  derivatives. 


4.  ABSORPTION 

Equation  (15)  and  Eqs.  (21)  and  (22)  give,  respectively,  the 
local  refractive  index  and  the  local  incident  power  intensity. 
Thus  the  amount  of  power  absorbed  per  unit  surface  area  is 
obtained  from  the  absorptivities  a,  and  ap,  respectively,  for 
the  s  and  p  waves  as  follows1: 

a,  ■  sin  20,,  sin  20in/sin2(0in  +  0r/), 


(23a) 


op  =  sin  20,,  sin  20in/sin2(0in  +  0,,)cos2(0in  -  0,,).  (23b) 

After  decomposing  the  different  trigonometric  functions  in 
Eqs.  (23)  and  substituting  (1  —  sin2  diJv2)1^2  and  sin  8\Jv, 
respectively,  for  cos  0,,  and  sin  0,,,  we  finally  get 

aa(Z,  <r,  «)  =  4  cos0in(i/2  -  sin20in)1/2/ 

[( i '2  ~  sin20in)1/2  +  cos  0in]2  (24) 


1/2 


A 


and 

ap(Z,  a,  w)  =  4x2  cos  0in( v2  -  sin20in)1/2/[(i<2  -  sin20in) 

+  cos0in]2[cos  0in(i<2  -  sin20in)1/2  +  sin20,„]w. 

(25) 

Equations  (24)  and  (25)  give  the  local  absorptivities  for  the  s 
and  p  waves  in  terms  of  Z  and  its  derivatives,  the  angles  of 
incidence  with  respect  to  a  reference  frame,  the  optical  aber¬ 
ration  u,  and  the  physical  properties  of  the  media.  Thus  the 
local  net  power  absorbed  is 


Ja.+P  =  +  apjp. 


(26) 


Integrating  Eq.  (26)  over  the  surface  area  exposed  to  the 
incident  rays  gives  the  total  absorbed  power. 


5.  SCATTERING 

Because  the  scattered  wave  front  propagates  in  a  nondissi- 
pative  medium,  interference  dominates  absorption.  There¬ 
fore  the  electric  and  magnetic  field  vectors,  and  not  the 
energy  fluences,  are  the  concerned  variables.  Equations 
(15)  and  (19)  give,  respectively,  the  local  refractive  index  and 
the  local  incident  field  vectors.  The  Fresnel  coefficients  of 
reflectance  R,  and  Rp,  respectively,  for  the  s  and  p  waves  are 
then  given  by1 

(27a) 


and 


R,  *  -sin(0in  -  0,,)/sin(0in  +  0,,) 


Rp  "  tan(0in  -  0,,)/tan(0in  +  0,,). 


(27b) 


In  a  way  similar  to  that  used  with  Eqs.  (24)  and  (25)  for  the 
absorptivities,  we  get  the  following  expressions  for  the  Fres¬ 
nel’s  reflection  coefficients: 

R,  =  -(O'2  -  sin20jn)1/2  -  cos  0,n]/[(*'2  -  sin20in)l/2  +  cos0in), 

(28) 

Rp  *  -R,/(cos  0in(i'2  •-  sin0in)1/2  -  sin20in).  (29) 


Equations  (28)  and  (29)  give  the  reflection  coefficients  for 
the  s  and  p  waves  in  terms  of  Z  and  its  derivatives,  the  angles 
of  incidence  with  respect  to  a  reference  frame,  the  optical 
aberration  u,  and  the  physical  properties  of  the  media. 

The  far-zone  scattered  electric  field  E,  at  a  point  in  the 
space  surrounding  the  irradiated  surface  is  obtained  by  the 
SCSI  as4 

E,(P)  =  -[ik,  exp(-i/i,fi)/4W?]Ap 

X  J[A  X  e, -  t )hp(h  X  h,)]  X  exp(ife,r  •  Ap)dS, 

(30) 

where  R  is  the  distance  from  the  origin  to  the  point  P,  hp  is 
the  unit  vector  in  the  direction  of  observation  and  dS  is  a 
surface  element.  The  surface  currents  in  Eq.  (30)  are  de¬ 
fined  by4 

A  X  e,  =  A  X  ej(l  +  R.)  +  ij( A  •  A,)(l  -  Rp)h‘p  (31) 


and 


A  X  h,  =  A  X  hj,(l  +  Rp)  ~  (A  •  A,)(l  -  RtWJn.  (32) 

All  terms  in  Eqs.  (31)  and  (32)  are  fully  determined  by  Eqs. 
(2),  (19),  (28),  and  (29).  Thus  the  scattered  field,  Eq.  (30), 
can  be  represented  explicitly  in  terms  of  Z  and  its  deriva¬ 
tives,,  the  angles  of  incidence  with  respect  to  a  reference 
frame,  the  optical  aberration  u,  and  the  physical  properties 
of  the  media. 

’ 

6.  APPLICATION  TO  A  SPHERICAL  WAVE 
FRONT  CONVERGING  ON  A  FINITELY! 
CONDUCTING  SURFACE 

A  spherical  wave  front  converging  on  a  conducting  surface 
haa  a  common  concern  in  many  optical  fields,  such  as  in 
inertial  confinement  fusion  and  atomic  cooling  by  focused 
laser  beams. 

A  spherical  wave  front  can  be  produced  by  an  aplanatic 
lens  as  follows.  Let  an  aberration-free  convex  lens  of  spheri¬ 
cal  surfaces  of  revolution  be  located  with  its  axis  of  revolu¬ 
tion  coinciding  with  the  z  axis.  A  plane  wave  propagating  in 
the  negative  z  direction  will  emerge  from  the  lens  as  a  spheri¬ 
cal  wave.  The  electromagnetic  field  vector  at  the  focal  vol¬ 
ume  of  the  lens  is  then  given  by,s 

E,  “  (£(/£,)  [  [  cosl/20  sin0A  exp(ifc,A,  •  r)d0dy>,  (33) 
Jo  Jo 

where  y  is  the  semi-aperture  angle  of  the  lens,  f  is  its  focal 
length,  and  Eo  is  the  electric  field  strength  at  the  lens.  The 
solution  of  Eq.  (33)  is  given  in  terms  of  the  diffraction  inte¬ 
grals  h,  and  I3  as15 


Exi  -  -iE0fk,{I0  +  /2  cos  2<p), 

(34a) 

Eyi  -  -iE0fk,I2  sin  2v>, 

(34b) 

Eti  =  -2 £„/*,/,  cos  ifi. 

(34c) 

Hxi  -  -iE0fklI2  sin  2v?/i), 

(34d) 

Hy ,  -  -iE0/A,(/0  -  /jcos  2 v>)/n, 

(34e) 
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and 

H,.  =  -2 E0fklIl  sin  v/h-  (34f) 

The  locally  polarized  waves  are  now  functions  of  the  diffrac¬ 
tion  integrals  and  are  obtained  by  substituting  the  E(  terms 
from  Eqs.  (34)  into  Eqs.  (19).  Equations  (21)  and  (22)  then 
give  the  local  energy  fluences  in  each  wave,  and  Eqs.  (24)  and 
(25)  give  the  absorptivities. 

Scattering  is  abo  treated  in  a  similar  way.  Equations  (31) 
and  (32)  give  the  surface  currents  in  terms  of  e,  and  h„  Eqs. 
(28)  and  (29)  give  the  reflection  coefficients,  and  Eq.  (30) 
gives  the  scattered  field. 

The  analogy  between  our  analysis  and  that  of  Ref.  2  is  as 
follows.  First,  the  far-field  amplitude  for  the  scattered 
field,  F'(6,,  <p,\6'it  *?,),  of  Ref.  2  can  be  expressed  from  our 
analysis  in  terms  of  the  surface  currents,  Eqs.  (31)  and  (32), 
as  will  be  shown  in  a  forthcoming  paper.  Second,  the  scat- 
terer  matrix  7V  of  Ref.  2,  which  contains  the  eight  surface 
integrals  of  spherical  vector  harmonics,  is  equivalent  to  the 
two  coordinate-dependent  reflection  coefficients  in  Eqs. 
(28)  and  (29),  but  the  inclusion  of  the  scatterer  shape  and 
physical  properties  in  the  R  terms  is  simpler  and  more  ex¬ 
plicit  than  that  in  the  7V. 


7.  CONCLUSION 

We  generalized  the  solution  for  a  plane  wave  incident  upon  a 
plane-conducting  surface10  to  the  case  of  arbitrary  surface 
and  wave-front  shapes.  Analytical  expressions  are  obtained 
for  the  local  refractive  index,  the  absorptivities,  and  the 
reflection  coefficients  for  a  finitely  conducting  arbitrary  sur¬ 
face.  The  explicit  inclusion  of  the  surface  function  and  its 
derivatives  into  our  solution  makes  it  possible  to  follow  the 
time  history  of  the  laser-solid  interaction  process  when  ther¬ 
modynamic  changes  take  place.  The  solution  obtained  is 
useful  in  computing  the  scattering,  as  well  as  the  extinction 
parameters,  for  a  wide  range  of  surface  geometries.  Our 
solution  is  also  simpler  than  the  methods  of  solution  by  the 
Laplace  transform,  the  Debye  potentials,  and  the  .vector 
wave  equation. 2.'6-17  Finally,  since  the  surface  integrals  of 
both  the  e  and  the  h  terms  are  devoid  of  spherical  vector 
harmonics,  the  computation  time  is  expected  to  be  greatly 
reduced  compared  with  that  of  the  method  described  in  Ref. 
2. 
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